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We investigate the excited states of the quasi-one-dimensional quantum antiferromagnets on hexagonal lat- 
tices, including the longitudinal modes based on the magnon-density waves. A model Hamiltonian with a 
uniaxial single-ion anisotropy is first studied by a linear spin- wave theory based on the one -boson method; the 
ground state thus obtained is employed for the study of the longitudinal modes. The full energy spectra of 
both the transverse modes (i.e., magnons) and the longitudinal modes are obtained as functions of the nearest- 
neighbor coupling and the anisotropy constants. We have found two longitudinal modes due to the non-collinear 
nature of the triangular antiferromagntic order, similar to that of the phenomenological field theory approach by 
Affleck. The excitation energy gaps due to the anisotropy and the energy gaps of the longitudinal modes without 
anisotropy are then investigated. We then compares our results for the energy gaps at the magnetic wavevectors 
with the experimental results for several antiferromagnetic compounds with both integer and non-integer spin 
quantum numbers. We also discuss the possible nonlinear effects (i.e., the so-called backflow corrections) in 
our microscopic approach. 

PACS numbers: 75.10.Jm, 75.30.DS, 75.50.Ee. 



I. INTRODUCTION 

The excitations of the quasi-one-dimensional (Id) Heisen- 
berg antiferromagnets systems have been studied extensively 
since Haldane predicted an energy gap in the excitation 
spectra of the isotropic integer-spin Heisenberg chains in 
1983 [1]. Now it is well established that there is an en- 
ergy gap separating the singlet ground state from the triplet 
lowest-energy-excitation states for the integer-spin Heisen- 
berg chains, contrast to the gapless excitation states of the 
half-odd-integer-spin Heisenberg systems [2, 3]. This theoret- 
ical prediction has been confirmed by Buyers et al [4] in the 
inelastic-neutron- scattering experiments on the quasi- Id anti- 
ferromangetic compound CsNiCl3. Some subsequent exper- 
imental investigations [4-9] and numerical calculations [10- 
14] also support Haldane 's prediction. 

At very low temperature, most of the quasi- Id antifer- 
romagnetic materials including CsNiCl3 show the three- 
dimensional nature with the classical magnetic order, and 
more interestingly, energy gaps at the magnetic wavevector 
have also been observed in many compounds [4]. For the case 
of CsNiCl3, the observed energy gap was initially explained 
by a uniaxial single-ion anisotropy but now it is widely ac- 
cepted that the gapped excited state belongs to one of the 
two longitudinal modes corresponding to the oscillations in 
the magnitude of the magnetic order of the quasi- Id hexago- 
nal systems, first proposed by Affleck based on a simplified 
version of Haldane 's theory [15, 16]. The gapped longitu- 
dinal modes are clearly beyond the conventional spin-wave 
theory which produces only the transverse excitations usually 
referred to as magnons. Later experimental study by Enderle 
et al. [17] using high-resolution polarized neutron scattering 
also confirms Affleck's proposal of the longitudinal modes, 
and contradicts to the spin-wave theory of two-magnon by 
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Ohyama and Shiba [18] or a modified spin- wave theory by 
Plumer and Caille [19]. There are also investigations of the 
longitudinal excitation states in other quasi- Id structures with 
the Neel-like long-ranged order at low temperature such as the 
tetragonal KCUF3 with s = 1/2 [20], where good agreements 
between the experiment and a theory based on a sine-Gordon 
field theory have been found for the energy gap at the mag- 
netic wavevector [21, 22]. More recently, a longitudinal mode 
was also observed in the dimerized antiferromagnetic com- 
pound TICUCI3 under pressure with a long-ranged Neel order 
[23]. 

We recently proposed a general microscopic many-body 
theory based on the magnon-density waves for the longitudi- 
nal excitations of spin-s quantum antiferromagnetic systems 
[24, 25]. In analogy to Feynmann's theory of the low-lying 
excited states in the helium-4 superfluid [26, 27], we iden- 
tify the longitudinal excitation states in a quantum antiferro- 
magnet with a Neel-like order as the collective modes of the 
magnon-density waves. In application to the quasi- Id tetrago- 
nal structure of KCuF 3 with s = 1/2, with no other fitting 
parameters than the nearest-neighbor coupling constants in 
the model Hamiltonian, we find our numerical results for the 
energy gap values at the magnetic wavevector are in general 
agreement with the experiments [28]. We hope that more ex- 
perimental results for the energy spectra at other wavevectors 
will be available for comparison. 

In this article, we extend our microscopic approach to 
the quasi- Id hexagonal quantum antiferromagnets such as 
CsNiCl 3 and RbNiCl 3 [17, 29, 30] both with spin-1 and 
CsMnI 3 with spin-5/2 [31]. The basal planes of these ma- 
terials are antieferromagnetic triangular lattice with the non- 
collinear magnetic order. Hence there are two possible lon- 
gitudinal modes in these hexagonal systems, rather than the 
single longitudinal mode of the bipartite systems such as the 
tetragonal KCUF3. Some preliminary results for the two di- 
mensional triangular model have been published [32]. We 
organize this article as follows. For completeness, we out- 
line the main results of the spin- wave theory for the quasi- Id 
model in Sec. 2, using the one-boson approach after two spin 
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rotations. We obtain the full spin- wave spectra as a function 
of the uniaxial single-ion anisotropy. To our knowledge, this 
anisotropy dependence of the spin- wave spectra has not been 
published before. We then apply our microscopic theory for 
the longitudinal excitations in Sec. 3, using the approximated 
ground state from the spin-wave theory. The energy gaps 
due to the anisotropy and the energy gaps of the longitudinal 
modes without anisotropy are then discussed in details. We 
compare our results for the energy gaps with the experimental 
results for the spin-1 compounds CsNiCl3 and RbNiCls and 
the spin-5/2 compound CsMnI 3 . We conclude this article by 
a discussion of the possible backflow corrections due to the 
nonlinear effects in the last section. 



spin rotations. Firstly, we rotate the local axes of all up-spins 
(blue) by 180° so that all spins along each chain align in the 
same down direction. This is equivalent to the transformation 



-ST 



Sf -> -Sf 



(2) 



for the first terms in Eq. (1), leaving the last two terms un- 
changed. The second rotation is on the hexagonal lattice of 




II. THE SPIN- WAVE THEORY OF THE ANISOTROPIC 
HEXAGONAL ANTIFERROMAGNETIC SYSTEMS 

The quasi- Id materials such as CsNiCl3 crystallize in the 
hexagonal ABX% structure with space group P63/mmc, 
where A is alkaline-metal cation, B is cation of the 3d group, 
and X is halogen anion. The magnetic ions B constructs the 
hexagonal lattice in the ab plane with adjacent spins forming 
angles of — 2tt/3, and antiparallel adjacent spins along the 
chain of the c axis as shown in Figs. 1(a) and (b). The lat- 

o 

tice constants of CsNiCl3 for example are a = 7.144 A and 

c = 5.90 A, and the magnetic moments are carried by Ni 2+ . 
The superexchange interaction between B (Ni 2+ ) ions is mod- 
eled by an N- spin Heisenberg Hamiltonian with a strong in- 
trachain interaction J and weak interchain interaction J' such 
as 



chain 



plane 

2J'^S;.S j + Z)^(Sf) 2 , (1) 



(id) 



where the notation indicates the nearest-neighbor cou- 
plings only and where we have also added an Ising-like single- 
ion anisotropy term with constant D(< 0). Most of the intra- 
chain couplings in ABX% compounds are antiferromagnetic 
such as CsNiCl3 or RbNiCl3 with easy single-site anisotropy, 
or CsMnBr3 and RbMnBr3 with hard anisotropy [33, 34]. 
These intrachain couplings can also be ferromagnetic (i.e., 
J < 0) as in CsNiF 3 [35, 36] or CsCuCl 3 [37]. We con- 
sider only the antiferromagnetic couplings here. Therefore, 
the classical ground state of each linear chain along the c axis 
(also denoted as y-axis) is a Neel state with alternating spin- 
up (blue) and spin-down (red) alignments as shown in Fig. 1 
(b). 

We consider a spin-wave theory for the Hamiltonian of 
Eq. (1) based on the one-boson approach by performing two 
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FIG. 1. The classical spin structure of the quasi- Id hexagonal an- 
tiferromagnets: (a) on the ab plane, and (b) the three-dimensional 
structure. 



the ab plane (or xz-plane) on the second and third terms of 
Eq. (1). Following Singh and Huse [38] and Miyake [39], for 
every triangle of the hexagonal lattices (see Fig 1(a)), we ro- 
tate the local axes of two spins along the classical direction in 
the xz-plane to align with that of the third spin. This is equiv- 
alent to the rotation of the i-sites of Eq. (1) by the following 
transformation 



Sf -> Sf, 

Sf ->Sf cos(0O 



■SfsmiPi), 
Sfsmft), 



(3) 



where 6i = Q z and Q z = (47r/3, 0, q z ) with Q z atq z = tt 
defined as the magnetic-ordering wavevector of the quasi- Id 
hexagonal systems. The Hamiltonian of Eq. (1) after these 
two transformations is given as 
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H = ~ o J E l S ? S te + S fSr +6 + 2SfSf +Q ] - - J' J2 [sw 



l,g 



l,Q> 



J l + Q' 



Si s l+e ,) 



-AS+S W + SrS+ +e ,) + 2 sin(6>, - l+e ,)(S;Sf +e , - SfSf +e ,)] + H D , 



(4) 



where I runs through all sites, g and q' are the nearest neigh- bor index vectors with coordination numbers z = 2 along the 
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chain and z' = 6 on the hexagonal basal planes respectively, 
and % D is the rotated anisotropy term. Care should be taken 
for the two rotations on this anisotropy term. The first rota- 
tion of Eq. (2) leaves it unchanged due to its quadratic form 
as mentioned before. In order to perform the second rotation 
of Eq. (3) involving rotations of the axes of the two spins to 
align with the axis of the third spin on the triangular planes, we 
rewrite the anisotropy term of Eq. (1) in the following equiva- 
lent, suitable form 



l,Q' 



3(^+0' 



(5) 



The transformation of Eq. (3) to the second term in Eq. (5) 
gives 



1,6' 

+ (Sl + Q' 



2 sin 2 6i +e > - cos sinfl| +e /(5f J _ n /5f_ 1 _ 



n+Q ,L >i+Q' 



Si+q>Si+6>) 



(6) 



We notice that this anisotropy form is different from the sim- 
ple form of Ref. [44] or that of Ref. [45] where the first term 
of Eq. (6) is missing. We believe that Eq. (6) is correct form 
suitable for the hexagonal systems. The energy gaps in the 
energy spectra due to this anisotropy term will be presented 
later. 

Using the canonical Holstein-Primakoff transformations, 
the spin operators are expressed in terms of a single set of 
boson operators a) and a as, 



(7) 



where / = y/l — a^a/2s a and s is the spin quantum num- 
ber. The Hamiltonian of Eq. (4) can then be written as, af- 
ter Fourier transformations of the boson operators with the 
Fourier component operators a q and a q and to the order of 

(25), 



H « Hq + Hi , 
where Hq is the classical energy, 



(8) 



H = -2JNs 2 +-3J'Ns 2 + -DNs 2 (l+2cos 2 <9+-sin 2 <9) 

3 s 

(9) 

with = 27r/3 and H 1 is given by the quadratic terms in the 
boson operators as 



with constants A q and B q defined by 



(10) 



A q = 4J + 6J'(1 



1 



>7d 



^D(l + 2cos 2 <9-sin 2 



B q = 4J cos q z + 9J'j q D sin 2 0, 



and 7^ defined as usual by 

1 iq-r 1 ( , Qx V3 

iq = —, e e = 3 { cos 2* + 2 cos y cos ~y% 

q' 

(12) 

The quadratic Hamiltonian Hi of Eq. (10) is diagonalized by 
the usual Bogoliubov transformation and can be written in 
terms of the new boson operators a q and a q as, 



(13) 



where AHq is the quantum correction to the classical ground 
state energy of Eq. (9), 

AH = -2 JNs - 3 J' TVs + -DNs{l + 2 cos 2 - sin 2 0), 

(14) 

and S q is the spin- wave excitation spectra, 



S q = sJA*-Bl 



(15) 




(11) 
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FIG. 2. (a) The first Brillouin zone of a quasi- Id hexagonal 
antiferromagnets. The points (0,0), (2tt/3, 2tt/\/3), (2tt/3,0), 
(4tt/3,0), (7r,7r/>/3), and (0,tt/>/3) all at ^ = 7r are denoted 
as Q / \K' \P' \Q,L' \O r respectively, and the similar points but at 
q z — are denoted as T,K,P,Q" \L,0 respectively, (b) The 
hexagonal Brillouin zone at q z —71 with some symmetry points in 
conventional notations for the quasi- Id systems. 

The first Brillouin zone of a quasi- Id antiferromagnet 
is ploted in Fig. 2, where the magnetic wavevector Q = 
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(47r/3,0,7r) is located at the corner of the hexagon and 
where other symmetry points in conventional notations are 
also illustrated. We plot the spin- wave spectra of Eq. (15) 
in Fig. 3 for CsNiCl3 using the experimental values J — 
0.345, J' = 0.0054 THz and negligible anisotropy D « 
[5, 15, 17, 42, 43]. We define the ratio of the two nearest- 
neighbor coupling constants as £ and for CsNiCl3, 

V 

i = — = 0.0157. (16) 

The spin- wave energy spectra with different polarizations are 
obtained by folding of the wavevectors. In Fig. 3, several 
branches along the symmetry direction of (0, 0, 77) , (77, 77, 1) 
and (1/3, 1/3, 1 + 77) are shown, where 77 is the reduced wave 
vector component in the reciprocal lattice unit (r.l.u) with 
q z = (2ttI/c) • (c/2) = irl, and 7 = l/3[cos7r/ + cos 2nk + 
cos27r(/i + k)]. Using Eq. (12) the moving in the param- 
agnetic Brillouin zone can be written as for q x = 4ttt] and 
Qz = + ^ and the corresponding symmetry directions to 
those in reciprocal lattice unit are (0, 0, tt + tttj), (47T77, 0, tt) 
and (47r/3, 0, tt + ^77) respectively. The three transverse spin- 
wave branches are obtained from Eq. (15) as follows. The 
y-mode has the polarization along the y-axis of the hexago- 
nal lattice where the quantum fluctuation is at q; the other two 
modes are found in the xz-plane by translating the wavevector 
by a magnetic wavevector as q — > (q ± Q) and are denoted as 
z£± respectively. 

As can be seen from Fig. 3, at the magnetic wavevector Q, 
the y-mode is gapless for zero anisotropy (D = 0). However, 
as mentioned earlier, an energy gap about 0.41 (2 J) has been 
observed by the neutron scattering experiments for CsNiCl3 
[4]. This energy gap can be reproduced in the y-mode exci- 
tation by introducing an anisotropy with D = —0.0285 using 
our approximation of Eq. (6), also plotted in Fig. 3. If we use 
the simple form from Ref . [44] corresponding to setting = 
in Eq.(6), the required anisotropy is reduced by a little more 
than half with the value D = —0.0141. Both of these values 
are now considered too large for CsNiCl3 which has negligi- 
ble anisotropy. The conclusion is that the observed gaps are 
not of the transverse spin-wave spectra, but belong to the lon- 
gitudinal modes, as first proposed by Affleck [15, 16]. 

Now we turn our attention to the order parameter. The long- 
range order of the quasi- Id hexagonal systems is given by the 
three sublattice-magnetizations with the same magnitude but 
different orientations as shown in Fig. 1, and it is clearly non- 
collinear, contrast to the collinear case of the bipartite sys- 
tems. In the spin- wave theory with one boson method as de- 
scribed above, the magnitude of the sublattice magnetization 
can be expressed as 



M 



(17) 



where the quantum correction p is the magnon density defined 
as the ground-state expectation value of the boson number op- 
erator 
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FIG. 3. The three spin- wave excitation spectra (in colors) for 
CsNiCl 3 with J = 0.345, J' = 0.0054 and D = THz, along the 
symmetry direction (0, 0, ty + tttj), (47t?7, 0, 7r) and 0, ty + tttj). 
Also included is the gapped ^/-mode (black, denoted as y') with 
D = —0.0285 using the anisotropy term of Eq. (6). The solid and 
dash with the blue color on the lines indicate the zx+-mode and zx~- 
mode respectively. 



with A q and B q defined by Eqs. (11). The numerical result 
of the magnon density for CsNiCl3 is p w 0.49 at D = 0, 
giving the sublattice magnetization M « 0.51. On the other 
hand, using slightly different parameter £ = 1.7xl0 -2 from 
Ref. [45], we obtain p = 0.48, giving M = 0.52. Both 
these results compare favorably with the experimental value 
of M = 0.53 [45]. As mentioned earlier, our microscopic 
analysis of the longitudinal modes is based on these mangon 
density fluctuations and there will be two such modes as dis- 
cussed in details in the following section. 



III. THE LONGITUDINAL MODES OF THE QUASI-1D 
HEXAGONAL ANTIFERROMAGNETS 



As mentioned before, the longitudinal excitations in a quan- 
tum antiferromagnetic system with a Neel-like long-ranged 
order correspond to the fluctuations in the order parameter. 
Using the fact that the quantum correction in the order pa- 
rameter is given by the magnon density p as discussed previ- 
ously in Eq. (17), the longitudinal modes can be considered 
as the magnon-density waves. By analogy to Feynman's the- 
ory on the low-lying excited states of the helium-4 superfluid 
[26], the longitudinal excitation states can be constructed by 
employing the magnon-density operators S z , in contrast the 
transverse spin- wave excitation states constructed by the spin- 
flip operators 5' ± [25]. The energy spectra of these longitudi- 
nal collective modes can then be easily derived by a formula 
first employed by Feynman for the famous phonon-roton spec- 
trum of the helium superfluid involving the structure factor of 
its ground state. 

More specifically, the longitudinal excitation state is ap- 
proximated by applying the magnon density fluctuation op- 
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erator X q to the ground state \^ g ) as, 



(19) 



where X q is given by the Fourier transformation of S z opera- 
tors, 



o, 



(20) 



with index I running over all lattice sites. The condition q > 
in Eq. (20) ensures the orthogonality to the ground state. The 
excitation energy spectrum in this linear approximation can be 
derived as [25], 



E(q) 



N(q) 



(21) 



where N(q) is given by the ground-state expectation value of 
a double commutator as 



N(q) = -([X- g ,[H,X q ]]) g , 



(22) 



and the state normalization integral S(q) is the structure factor 
of the lattice model 

S(q) = (X_ q X q ) g = ±J2 ei<l - (r, ~ r "HSiS[,) g . (23) 



The notation (. . . ) g in Eqs. (22) and (23) indicates the ground- 
state expectation. We have applied these formulas to the bipar- 
tite quasi- Id antiferromagnetic systems such as KCUF3 [28]. 
For the hexagonal lattice systems as discussed here, we ex- 
pect that there are two longitudinal modes due to the non- 
collinear nature of the order parameter on the triangular basal 
plane. Within the one-boson approach after the two spin ro- 
tations as employed here, the two longitudinal modes with 
xz-polarizations of the hexagonal systems can be obtained by 
folding of the wavevectors in the energy spectra of Eq. (21), 
in similar fashion to the one-boson spin-wave theory as dis- 
cussed in Sec. II and also to that of Ref. [15]. 

Using the Hamiltonian of Eqs. (4), it is straightforward to 
derive the following double commutator with zero anisotropy 
(i.e. D = 0) as 

N(q) = 2 Js ^(1 + cos q z )g e + \j' [3(1 + lq)~9 e > 

Q Q' 

~ (1 -Jq)g f e>\i 

(24) 

where j q is as defined in Eq. (12) and the transverse correla- 
tion functions g r and g' r are defined respectively as 



2s 



all independence of index I due to the lattice translational 
symmetry. Their Fourier transformations are obtained as, us- 
ing the approximation of the linear spin- wave theory for the 



ground state in the one-boson approach discussed in Sec. 2, 



B„ 



9q= ~ 



2 Jai-bi 




1 , (26) 



with A q and B q as given before by Eqs. (11). The numerical 
results are g e = 0.669, g Q , = 0.056 and g' Q , = -0.0076. As 
can be seen, N(q) is dominated by g e . The structure factor 
within the linear spin-wave approximation is given by 



S(q) = P + jf^Pq'Pq+tf + ^ J2^' 



N ■ 



9q+q 



(27) 



where p — ^ J2 q p q is the magnon density of Eq. (18), with 
definition 




A\-Bl 



(28) 



We first discuss the general behaviors of the longitudinal 
spectrum of Eq. (21) as a function of the ratio of the two 
nearest-neighbor coupling constants, £ = J' / J . In the limit 
£ — >> 0, the Hamiltonian of Eq. (1) becomes the pure Id sys- 
tems; the longitudinal spectrum is gapless and identical to 
the doublet spin-wave spectra thus forming a triplet excita- 
tion state as discussed in details Ref. [28]. This demonstrates 
the limitation by the spin- wave ground- state employed, par- 
ticularly when applied to the integer-spin Heisenberg chain 
where the Haldane gap is expected as discussed in Sec. I. In 
the other limit, £ — » oo, the Hamiltonian is a pure triangular 
antiferromagnet with the quasi-gapped longitudinal modes as 
discussed in details in our previous paper [32], similar to the 
case of the square lattice model. 

For the quasi- Id materials with intermediate values of £, 
the spin- wave ground state is a reasonable approximation. 
We obtain nonzero energy gaps for the longitudinal excitation 
spectra of Eq. (21). As discussed before, following Affleck 
[15, 16], two longitudinal modes for the quasi- Id hexagonal 
antiferromagnets can be obtained by folding of the wavevec- 
tor. We denote one as L_ with the spectrum E(q — Q) and the 
other as L + with the spectrum E(q + Q). We plot these two 
longitudinal spectra together with the three spin- wave spectra 
of Eq. (15) in Fig. 4 near the magnetic wavector Q for the 
compound CsNiCl3. Our numerical result for the energy gap 
of the lower longitudinal mode L_ at Q is (0.96) 2 J, more 
than twice of the experimental results of 0.41 (2 J). We also 
notice that the upper mode L + is higher than the L_ mode 
by about (0.18)2 J at Q. As will be discussed in more details 
later, we believe that the higher values of our energy gap are 
due to the linear approximation in our analysis. We also plot 
the L_ mode along the the path Q' K' PQLO of the hexag- 
onal Brillouin zone in Fig. 5 together with the spin- wave y 
and zx- modes. As can be seen, over the whole spectrum, 
the longitudinal mode in our linear approximation does not 
change much. 

For the compound RbNiCl3 also with 5 = 1, using the ex- 
change parameters J = 0.485 and J' = 0.0143 THz with a 
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ber s = 5/2 and the nearest-neighbor coupling constants 
J = 0.198 and J' = 0.001 THz and negligible anisotropy 
[31]. This system is very close to the pure Id system with 
a very small ratio £ = J' / J w 0.005. The linear spin- wave 
theory may be a poor approximation for such a system. Never- 
theless, with the similar analysis as before based on the spin- 
wave ground state, we obtain a gap value of 0.64 THz for the 
lower L_ mode at the magnetic wavevector Q. This is much 
larger than the experimental value of about 0.1 THz by Har- 
rison et al [31], which was used to fit a modified spin- wave 
theory by Plumer and Caile [19]. Clearly, for such systems as 
CsMnl3, we need a better ground state than that of the spin- 
wave theory in our analysis, in additional to the nonlinear ef- 
fects in the excitation state operator as mentioned before for 
CsNiCl3 and RbNiCl3 compounds. We will discuss further in 
the next section. 



FIG. 4. The longitudinal modes L± as derived from Eq. (21) together 
with the spin- wave y- and zx± modes as derived from Eq. (15) for 
CsNiCb along the symmetry direction (0,0, ty + irrj), (47r?7,0,7r) 
and (4^, 0, 7r + tttj). 
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FIG. 5. The longitudinal mode L_ along the path Q' K' PQLO of 
the hexagonal Brillouin zone of Fig. 2(b) together with the spin-wave 
y and zx~ modes for CsNiCl3. 



larger ratio £ = J' / J = 0.0295 [49], we obtain similar lon- 
gitudinal modes as those of CsNiCl3. The numerical result 
for the energy gap of the L_ mode is 1.16 THz at the mag- 
netic wavevector, again about twice of the experimental value 
of 0.51 THz. We like to point out that there is some difficulty 
in fitting of Affleck's model with the experimental results for 
RbNiCl 3 [16,49]. 

Finally we turn to the longitudinal modes for the non- 
integer-spin quasi- Id hexagonal systems. The superexchange 
interactions in the hexagonal compound CsMnl3 can be de- 
scribed by the Hamiltonian of Eq. (1) with spin quantum num- 



IV. DISCUSSION 

In this paper, we have investigated the excitation states of 
the quasi- Id hexagonal systems as modeled by the anisotropic 
Heisenberg Hamiltonian with only the nearest-neighbor cou- 
plings. We have obtained the three spin- wave modes and two 
longitudinal modes. The energy gaps due to the anisotropy 
and the energy gaps of the longitudinal modes at the mag- 
netic wavevector are investigated and compared with the ex- 
perimental results for several quasi- Id hexagonal compounds. 
We like to emphasize that our analysis applies to both inte- 
ger and non-integer spin systems and there are no other fitting 
parameters than the nearest-neighbor coupling constants and 
the anisotropy parameter in the model Hamiltonian. For the 
compounds CsNiC^ and RbNiCl3, our estimate for the en- 
ergy gaps at the magnetic wavevector is about twice of the 
experimental values; for the compound CsMnl3 which is very 
close to the pure Id model, our estimate is much larger than 
the experimental value. 

Our larger energy gap values than the experimental results 
are perhaps due to the linear approximation in our excita- 
tion operator X q of Eq. (20), where the couplings between 
magnons and the longitudinal modes have been ignored. It is 
interesting to note that for the case of the phonon-roton spec- 
trum of helium superfluid, a similar linear theory produces an 
energy gap near the roton wavevector with a value of about 
twice of the experimental results, and better agreement with 
the experimental results is obtained only after inclusion of 
the nonlinear effects (or the so-called backflow corrections) 
[27, 51]. We therefore believe that such nonlinear effects may 
also be significant in the quasi- Id hexagonal quantum ant- 
ferromagnets and deserve further investigation. Furthermore, 
particularly for the compound CsMnl3 where the interchain 
coupling is particularly weak, a better ground state than that 
of the linear spin- wave theory is needed. A more sophisticated 
many -body theory such as the coupled-cluster method, partic- 
ularly its recent variational version [52, 53], may provide such 
improvement. 
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